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The Kondo Lattice Model in Infinite Dimensions II. 
Static Susceptibilities and Phase Diagram 
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Magnetic and charge susceptibilities in the Kondo lattice are derived by the continuous-time 
quantum Monte Carlo (CT-QMC) method combined with the dynamical mean-field theory. 
For a weak exchange coupling J and near half filling of the conduction band, antiferromag- 
netic transition occurs as signalled by divergence of the staggered magnetic susceptibility with 
lowering temperature. With increasing J, the Kondo effect suppresses the divergence, and the 
critical value of J agrees well with Doniach's estimate which considers the RKKY interaction 
as competing with the Kondo effect. For low density of conduction electrons, a ferromagnetic 
ordering is observed where Doniach's estimate does not work. Around quarter filling, a charge- 
density-wave (CDW) transition is found. The CDW is interpreted from the strong-coupling 
limit in terms of effective repulsion between Kondo singlets. 

KEYWORDS: Kondo lattice model, charge-density wave (CDW), Doniach's phase diagram, continuous- 
time quantum Monte Carlo method (CT-QMC), dynamical mean-field theory (DMFT), 
hypercubic lattice 



1. Introduction 

Due to the strong local correlations, 4/ electrons of 
rare-earth ions exhibit itinerant as well as localized prop- 
erties. A single magnetic impurity in a metal gives rise 
to the Kondo effect, and the RKKY interaction works 
between a couple of moments. Consequently, rare-earth 
materials, considered as systems with periodic array of 
the localized moments, show a variety of phases. The sim- 
plest model to describe the above situation is the Kondo 
lattice model given by 

c^ + jysi-^, (i) 
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where Si denotes the localized spin at i site, and cr\ is 
defined by ct\ = J2aa' c i a <T ^' c icr' ■ With the spectrum e*. 
fixed, the system described in eq. (1) is specified by two 
parameters: a number of conduction electrons per site, 
n c , and the coupling constant J. 

As for the one-dimensional system, the ground-state 
properties of the Kondo lattice have been revealed based 
on the diagonalization and density matrix renormaliza- 
tion group methods. 1,2 Although strong quantum fluctu- 
ations suppress an antiferromagnetism in one dimension, 
a ferromagnetism can survive since the total spin com- 
mutes with the Hamiltonian. Actually, the ferromagnetic 
phase is stabilized in the low-carrier and strong-coupling 
regime, while the other region is paramagnetic. 

Although an overall picture has been clarified in one 
dimension, these results cannot remain valid in higher 
dimensions, where antiferromagnetism takes place. Vari- 
ational methods have been used for a ground-state phase 
diagram that includes ferromagnetic and antiferromag- 
netic phases at weak coupling. 3,4 More reliable results 
have been obtained in two dimensions with use of a quan- 
tum Monte Carlo simulation. However, the simulation is 
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restricted to the half-filled case of the conduction band 
because of the notorious minus sign problem. 5 On the 
other hand, higher dimensions can be treated by means 
of the dynamical mean-field theory (DMFT). With use of 
the numerical renormalization-group method, the single- 
particle spectrum 6 and the magnetic phase 7 have been 
derived. A fundamental problem is how the Kondo ef- 
fect works in the presence of many spins, which number 
is comparable to or even larger than that of conduction 
electrons. 8 

In a previous paper, 9 we have derived formulae to eval- 
uate spatial dependences of two-particle correlations in 
the framework of the DMFT. In the DMFT, the strong 
local correlations are rigorously taken into account. 10 
Therefore, we can discuss the magnetic and charge insta- 
bilities against formation of the coherent quasi-particles. 
We employ the recently developed CT-QMC algorithm 
to solve the effective impurity Kondo model. 11-13 With 
these framework, numerically reliable solutions can be 
obtained. We will demonstrate existence of antiferro- 
magnetic and ferromagnetic phases as well as the CDW 
phase, each of which is stabilized depending on the 
conduction-electron density n c . 

This paper is organized as follows. In the next sec- 
tion, we review the properties of the tight-binding band 
in the infinite-dimensional hypercubic lattice. Numeri- 
cal results of the susceptibilities are given in §3-6 with 
varying the occupation number of the conduction band. 
From these results, a ground-state phase diagram is con- 
structed in §7. Discussions will be finally given in §8. 

2. Infinite-Dimensional Hypercubic Lattice 

The cubic lattice can be extended to higher dimen- 
sions, which is referred to as the hypercubic lattice. 
The tight-binding band with a nearest-neighbor hopping 
has the nesting property at half-filling with the vector 



2 J. Phys. Soc. Jpn. 



Full Paper 



J. Otsuki et al. 



Q = (n, ■ ■ ■ , 7r). The nesting may give rise to instabilities 
of the Fermi surface resulting in, for example, magnetic 
order or CDW. In this section, we give a brief summary of 
properties of the infinite-dimensional hypercubic lattice, 
following rcfs. 10 and 14. 

We consider a hypercubic lattice in dimension d with 
the nearest-neighbor hopping t. We scale the hopping as 
t = D/2\/2d to be applicable to d = oo. The energy 
dispersion is then given by 



D 



(k = - 



2d z 



^ cos ki 



(2) 
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where we take the lattice constant unity. If the vector 
k satisfies = 0, which is a condition of the half-filled 
Fermi surface, then the vector k + Q also gives ek+Q = 0. 
Consequently, in any dimension, the hypercubic lattice 
exhibits the complete nesting property at half filling. 

The density of states p(e) = N^ 1 J2 k S(e — e&) with 
iV being the number of lattice sites can be evalu- 
ated by means of the Fourier transform. In the infinite- 
dimensional limit, d = oo, p(e) becomes the Gaussian 



(3) 



Therefore, the energy band has no cutoff in infinite di- 
mensions. We choose the unit D = 1 in this paper. 

To see characteristics of responses in the hypercubic 
lattice, we evaluate the polarization function of free con- 
duction electrons, 11°, defined by 



n" 



1 \ ~> fk+q — fk 

N ^ e k - e k + q ' 
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where fk = 1/(1 + c^ h ^^/ T ). The staggered component 
at half filling diverges logarithmically with decreasing 
temperature due to the complete nesting, e& = £k+Q- 
In the limit of d = oo, the q-dependence appears only 
through the quantity 77(g) defined by 14 



1 d 
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Figure 1 shows filling dependences of 11° at T = 0. At 
n c = 0.8, a distinct maximum is found close to q = Q 
(77 = —1). At a density less than n c — 0.6, on the other 
hand, the uniform component (77 = +1) is favored rather 
than the staggered one (77 = —1). 

Even taking the interaction into account, the q- 
dependence of the two-particle correlation is represented 
by 77(g) as long as the self-energy is local. Especially for 
77 = +1 and —1, the summation over k can be performed 
analytically. In this paper, we compare the uniform and 
staggered components without considering possibility of 
incommensurate ordering. 

3. Properties at Half Filling 

We now present numerical results computed in the 
DMFT and the CT-QMC. Let us begin with the half- 
filled band. Figure 2(a) shows temperature dependences 
of the static susceptibilities of localized spins, xLc X q =o 
and Xq, for J = 0.4. For comparison, we also plot the 




(q = 0) 



(q = Q) 



Fig. 1. The polarization function 11° of free conduction electrons 
at T = in the d = 00 hypercubic lattice. 



susceptibility in the single-impurity system, Xi™ P ' ) - At 
high temperatures, the susceptibilities simply exhibit lo- 
cal fluctuations following the Curie law, x q — C/T with 

C = 1/4. The local component x{ oc takes a maximum 
at about T = 0.15. This temperature dependence is con- 
trastive to Xi™ P ' ) 1 which monotonously increases with de- 
creasing temperature. The reduction in the lattice sys- 
tem is due to a formation of an energy gap as shown 
later. Compared to the local susceptibility, the staggered 
component is largely enhanced while the uniform one 
is suppressed. The large antiferromagnetic fluctuation is 
ascribed to the nesting of the Fermi surface. 

We show, in Fig. 2(b), the magnetic and charge suscep- 
tibilities of conduction electrons, x c ' 7 (7 = ma g>chg)- At 
high temperatures, the susceptibilities are approximately 
given by x°' 7 — 2C 7 IT ( , where C 7 corresponds to C for 
the magnetic channel and unity for the charge channel. 
The temperature dependence of magnetic susceptibilities 
is similar to that of the localized spins for each case of 
X c ' mag , while the charge susceptibilities decrease at low 
temperatures. In particular, the uniform charge fluctua- 
tion Xq=o S eventually vanishes. This temperature depen- 
dence indicates formation of an energy gap, namely the 
Kondo insulator. 

For smaller coupling, the characteristic energy be- 
comes smaller with enhanced antiferromagnetic fluctu- 
ation. The Kondo ground state then becomes unstable 
against the RKKY interaction. Figure 3(a) shows tem- 
perature dependence of the inverse susceptibilities of the 
local spins for J = 0.2 at T < 0.2. We can recognize 
a divergence of the staggered magnetic susceptibility at 
Taf — 0.03, indicating the antiferromagnetic transition. 
Since these calculations assume the paramagnetic effec- 
tive medium at all temperatures, the negative suscepti- 
bility has been obtained below the transition tempera- 
ture. In the present approximation based on the DMFT, 
the singularity in the susceptibility has no influence on 
the self-consistency of the effective medium. 

Near the transition temperature, the antiferro- 
magnetic susceptibility follows the Curie-Weiss law, 
Xq(T)/C oc (T - Taf) -1 - This behavior originates in 
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Fig. 2. Temperature dependence of (a) the local (x{ oc )> uniform 
(Xq=o)> ana - staggered (Xq) components of static susceptibili- 
ties of local spins, and (b) the corresponding susceptibilities of 
conduction electrons. The parameters are chosen as J = 0.4 and 
n c = 1. 



the mean-field treatment of the intersite interaction in 
the DMFT, which is justified in infinite dimensions. The 
Curie- Weiss law has been reported also in the infinite- 
dimensional Hubbard model. 15 The proportionality fac- 
tor in the Curie- Weiss law is about 0.7 for J = 0.2, which 
demonstrates a reduction of the local moment due to 
the Kondo effect. Figure 3(b) shows the corresponding 
conduction-electron susceptibilities. The staggered com- 
ponent of the magnetic channel diverges at Taf affected 
by the corresponding fluctuation of the local spins. The 
charge fluctuations remain finite in contrast with the 
Kondo insulating state. 

To see the dependence of the transition temperature 
on J, we turn our attention to the staggered component 
Xq- F° r small- J case such that J < 0.2 in Fig. 4(a), 
the transition temperature Taf monotonously increases 
with J. The increased slope of 1 /Xq indicates reduction 
of the moment. The reduction factor for the Curie con- 
stant changes from about 0.9 at J = 0.05 to about 0.7 
at J = 0.20. As J becomes larger than 0.2 in Fig. 4(b), 
the susceptibility deviates from the Curie- Weiss law and 
tends to saturate. Following the saturation of \q, the 
transition temperature turns to decrease, and at J = 0.27 
we cannot see the divergence at finite temperatures. 



Fig. 3. Temperature dependence of inverse static susceptibilities 
for J = 0.2. 



We summarize the foregoing results on the antifcr- 
romagnctic instabilities at half filling. Figure 5 shows 
the antiferromagnetic transition temperature Taf as a 
function of J. For comparison, we show the transition 

f MF^ 

temperature T AF evaluated in the mean-field theory, 
cq. (A-4). At weak coupling of J = 0.05, the mean-field 
theory turns out to reproduce the DMFT results quanti- 
tatively. This follows from the fact that the DMFT deals 
with the intersite interactions as the mean field. As J 
increases, the transition temperature in the DMFT falls 
below the mean-field results owing to the development of 
the local correlation, which screens the local spins from 
the surrounding spins. The maximum of Taf is found to 
be about 0.03 at J ~ 0.22. After taking the maximum, 
Taf rapidly decreases around J ~ 0.25 and eventually 
vanishes at J c ~ 0.27. 

Following Doniach, 16 we estimate the critical value 
of the coupling constant by comparing the Kondo tem- 
perature Tk and the RKKY interaction. We invoke the 
single- impurity system to evaluate Tk- From the static 
susceptibility ^( im P) at zero temperature, Tk is defined 
by T K imp) = C/x (imp) - We also estimate the Kondo tem- 
perature with use of the following analytic expression 
derived by means of the scaling theory: 17 

^scaling) = £ v /2J-J7j7jex P r ' 



2Jp(p) 



(6) 



4 J. Phys. Soc. Jpn. 



Full Paper 



J. Otsuki et al. 



*5? 




Fig. 4. Temperature dependences of the inverse staggered suscep- 
tibility 1/Xq °f localized spins for several values of J. 



where \i denotes the chemical potential for J = 0. While 
this expression is originally derived for the constant den- 
sity of states, eq. (6) holds well also for the Gaussian den- 
sity of states as long as Tj< is much smaller than the band 
width. The Kondo temperature is plotted in Fig. 5 to- 
gether with the mean-field transition temperature T^f F ^ 
evaluated from eq. (A-4). It turns out that T^p 1 ^ crosses 
T£ mp) and T^ caling) between J = 0.25 and J = 0.30. 
Hence, comparison of the energy scales reproduces fairly 
well the critical coupling constant. These numerical re- 
sults support Doniach's picture for competition between 
the Kondo singlet and the magnetic ordering due to the 
RKKY interaction. 

As shown in Fig. 4(b), the staggered susceptibility 
around the critical coupling J c saturates at low temper- 
atures. Hence, the expansion of Xq(^) 1 with respect to 
T does not include a linear term near J c . We approximate 
this behavior as 



X q(T, jy 1 =aT 2 + b(J-J c ). 



(7) 



Equating Xq{Taf,J) 1 to 0, we obtain the transition 
temperature near J c as 

T AF (J)cx(J c - J) 1 ' 2 . (8) 

Consequently, the slope of Taf(J) diverges at J c . This 
feature can be seen in Fig. 5. 
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Fig. 5. Antiferromagnetic transition temperature Taf as a func- 
tion of J at half filling. T^ F ^ is the transition temperature 
evaluated in the mean-field theory. 7^ lmp ) and j^ scalm s) arc tnc 
Kondo temperature determined from the low-temperature static 
susceptibility in the single-impurity system and from the formula 
in the scaling theory, eq. (6), respectively. 




Fig. 6. Temperature dependence of the inverse staggered suscep- 
tibility 1/Xq °f localized spins near half filling with J = 0.2. 



4. Doping on Insulating States 

We have so far examined the half filling case n c = 1, 
and have observed two insulators: the antiferromagnetic 
insulator and the Kondo insulator. Let us next exam- 
ine the filling-dependence of these two phases. We ex- 
pect reduction of the antiferromagnetic fluctuation, and 
crossover to the metallic state from the Kondo insula- 
tor. Because of the particle-hole symmetry of the Kondo 
coupling, n c and 2 — n c are equivalent. Henceforth, we 
investigate hole-doping n c < 1 only. 

We first examine the antiferromagnetic phase. Fig- 
ure 6 shows temperature dependences of the inverse of 
the staggered susceptibility at J = 0.2 for several val- 
ues of n c near half filling. Doping on the half-filled band 
suppresses the antiferromagnetic fluctuation. This is due 
to the change of the nesting condition in the Fermi sur- 
face. At n c = 0.88, the antiferromagnetic transition dis- 
appears, and the paramagnetic ground state is stabilized. 
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Fig. 7. Antifcrromagnctic transition temperature Taf as a func- 
tion of J for several values of n c . 
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Fig. 8. Temperature dependence of the uniform charge suscepti- 
bility Xq'=o S °f conduction electrons near half filling for J = 0.4. 



In Fig. 7 we show an antiferromagnetic transition tem- 
perature near half filling. The phase boundaries seem to 
be almost scaled by J. At n c = 0.80, we have observed 
no transition down to T = 0.005. 

We turn our attention to the Kondo-insulator. Figure 8 
shows filling dependences of the uniform charge suscep- 
tibility, Xq=(^> f° r J = 0-4- The charge susceptibility is 
sensitive to the doping, which gives rise to particle-hole 
excitations and a finite charge susceptibility at absolute 
zero. Similar non-monotonous behavior of Xq=o S (^) nas 
been observed in the one-dimensional Kondo lattice. 18 
These two energy scales also appear in the single-particle 
excitation spectrum, which shows an energy gap of order 
Tk and the coherent state at low energies. 19 

5. Ferromagnetism at Low Carrier Densities 

At low density, the conduction electrons do not have 
the nesting property. The polarization function 11® in 
Fig. 1 exhibits a flat structure at wide range of q. 
The maximum occurs at q — rather than Q. The 
RKKY interaction accordingly favors ferromagnetic or- 
dering rather than antiferromagnetism at low carrier 




0.05 0.1 0.15 0.2 



T 

Fig. 9. Temperature dependence of inverses static susceptibilities 
for J = 0.4 and n c = 0.2. 

density. 

Let us compare the local susceptibilities in the single- 
impurity and periodic systems. If the number of conduc- 
tion electrons are less than the number of local spins, it 
is impossible to screen all the spins even in the strong 
coupling limit J = oo. This is in strong contrast with 
the impurity system. Figure 9(a) shows temperature de- 
pendence of inverse magnetic susceptibilities for J = 0.4 
and n c = 0.2. The inset shows that the susceptibility in 
the periodic system continues to grow in the temperature 
range where the impurity susceptibility is saturated. It 
turns out from Fig. 9(a) that the susceptibility of the 
Kondo lattice follows the Curie- Weiss law in a wide tem- 
perature range, which indicates an existence of the al- 
most free moment. In contrast, the impurity model has 
a large deviation from the Curie- Weiss law due to the 
Kondo screening. 

As expected from the behavior of I±° in Fig. 1, the 
uniform component is larger than the staggered one. The 
uniform susceptibility eventually diverges at a finite tem- 
perature, indicating the ferromagnetic transition. The 
corresponding susceptibility x c of conduction electrons 
is shown in Fig. 9(b). The ferromagnetic component di- 
verges following the behavior of local spins. 

Let us examine the filling-dependence of the ferromag- 
netic transition temperature at low carrier densities. Fig- 
ure 10 shows the inverse uniform susceptibility for vari- 
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Fig. 10. Temperature dependence of inverse uniform susceptibil- 
ity C/Xq =0 for J = 0.4 and several values of n c . 



ous values of n c up to 0.68 with J = 0.4. We can clearly 
see suppression of the ferromagnetic correlation with in- 
creasing n c . We estimate the ferromagnetic transition 
temperatures Tf by extrapolating the data at T — 0.005 
and 0.010 according to the Curie- Weiss law. At the low- 
carrier-density regime, such an extrapolation is expected 
to be applicable reasonably. Figure 11(a) shows T-p as a 
function of n c for J = 0.3, 0.4 and 0.5. As n c increases, Tp 
decreases and turn to negative values. For J = 0.4, the 
ferromagnetic transition disappears at about n c = 0.64. 

Figure 11(b) shows the transition temperature T F as 
a function of J at n c = 0.2, 0.4 and 0.6. After taking a 
maximum against J, Tp decreases gradually, in contrast 
with the antiferromagnetic case shown in Figs. 5 and 7. 
At n c = 0.2, the fcrromagnetism persists beyond J = 
1.3. If we compare Ik in eq. (6) and T^ MF ) in eq. (A-4), 
we obtain the critical coupling constant as J ~ 0.6 for 
n c = 0.2. Hence, Doniach's picture does not hold in the 
low-carrier-density regime, where the carrier number is 
insufficient to screen all local spins. 

6. CDW at Quarter Filling 

In this section, we show results for quarter filling 
n c = 0.5. Figure 12(a) shows inverse magnetic suscep- 
tibilities of localized spins at J = 0.3. The local spins 
have the ferromagnetic correlation. An extrapolation re- 
veals that the divergence of the ferromagnetic fluctuation 
takes place at the temperature of order 10~ 3 . 

We show the corresponding conduction-electron sus- 
ceptibilities in Fig. 12(b). The staggered component of 
the charge susceptibility Xq c1is diverges at about T ~ 
0.014, indicating an ordering into the chargc-density- 
wave (CDW) state. Concerning the magnetic channel, 
the uniform component is favored following the corre- 
sponding fluctuation of the local spins. For stronger cou- 
pling constant, J = 0.6, The ferromagnetic correlation 
remains toward divergence (Fig. 12(c)). On the other 
hand, the increase of the CDW fluctuation is suppressed 
at low temperatures as shown in Fig. 12(d). 

We concentrate our attention on the CDW instabil- 
ity. We first vary J keeping n c . Figure 13 shows l/Xq Chs 
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Fig. 11. Ferromagnetic transition temperature Tp (a) as a func- 
tion of n c for J = 0.3, 0.4 and 0.5, and (b) as a function J for 
n c = 0.2, 0.4 and 0.6. 



as a function of temperature in logarithmic scale. As J 
decreases, the transition temperature decreases, and for 
J < 0.18 no divergence is found down to T = 0.005. 
However, we cannot exclude the CDW ordering at these 
parameters, since the energy scale seems to decrease ex- 
ponentially against J as shown in Fig. 13(a). On the 
other hand, a boundary on the strong-coupling side ex- 
ists as is clear from Fig. 13(b). The susceptibility tends 
to saturate with increasing J. The CDW transition tem- 
perature Tcdw is plotted in Fig. 14. The maximum of 
Tcdw is about 0.016 at J - 0.36. The CDW transition 
is completely suppressed at J > 0.47. 

We next shift n c from quarter filling. Figure 15 shows 
i/Xq C1iS at J = 0.3 for various fillings around n c = 0.5. 
The transition temperature decreases on both sides off 
quarter filling. This result seems to be reasonable since 
the quarter filling is favorable for the CDW state at 
q = Q, where the conduction electrons tend to align at 
alternate sites. The mechanism of the CDW instability 
will be discussed in detail in §8 based on the strong- 
coupling expansion. 

7. Ground-State Phase Diagram 

Figure 16 shows a ground-state phase diagram against 
J / D and n c estimated from the finite-temperature calcu- 
lation. The boundaries arc determined through scanning 
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Fig. 12. Temperature dependence of the inverse static susceptibilities at quarter filling for J = 0.3 in (a) and (b), and J = 0.6 in (c) 
and (d). 



the inverse susceptibilities with varying either J or n c , 
as in Figs. 4, 10, 13 and 15. The error bars are parallel 
to the axis of the scanned parameter. 

Near half filling, antiferromagnetic and paramagnetic 
phases divide into weak- and strong-coupling regimes. 
The boundary agrees roughly with the condition of com- 
parable energy scales of the RKKY interaction and the 
Kondo effect, which demonstrates the validity of Do- 
niach's picture. The paramagnetic state at half filling is 
the Kondo insulator, as revealed by the uniform charge 
susceptibility in Fig. 2. 

The ferromagnetic phase appears for n c < 0.6, which 
is naturally understood from the RKKY interaction. The 
polarization function takes a lager value at q = rather 
than q — Q as shown in Fig. 1 . The ferromagnetic phase 
is enlarged as compared to the mean-field estimation of 
Doniach. The ferromagnetic phase is stable beyond J = 
1.3 at n c = 0.2. 

The CDW state at q = Q arises around the quar- 
ter filling n c = 0.5 with a range of about ±0.06. We 
could not accurately determine the boundary on the 
weak-coupling side because of the exponentially decreas- 
ing transition temperature. Hence, possible boundary is 
represented by the dashed line in Fig. 16. Since the CDW 
state still keeps a part of the entropy of localized spins, 
the magnetic ordering is expected at temperatures lower 
than the CDW transition. The magnetic susceptibilities 



at temperatures higher than the CDW ordering show a 
tendency toward the ferromagnetism. In fact, the CDW 
phase is surrounded by the ferromagnetic phase. There- 
fore, it is likely that the CDW state coexists with the fer- 
romagnetism in the ground state. For couplings smaller 
than J ~ 0.14, the energy scale of the charge fluctuation 
is smaller than that of the ferromagnetism. Hence, the 
CDW ordering is unlikely to occur in the small- J case. 

Around n c = 0.8, no divergence of the susceptibilities 
has been observed for any J down to 0.1. In view of the 
fact that the Kondo temperature is exponentially small 
in the weak coupling regime, it is highly probable that an 
instability actually takes places at certain ordering vector 
other than q = and Q. The polarization function TI^ of 
free electrons in fact has a maximum between q = and 
Q as shown in Fig. 1. For n c — 0.6 and 0.8, the peak is 
located close to q = Q. Hence, we reasonably expect the 
incommensurate order. The ordering vector q can change 
successively, so that the boundaries will be complicated 
between the ferromagnetic and incommensurate as well 
as between the incommensurate and antiferromagnetic 
phases. 

8. Summary and Discussion 

Insulating phases at half filling 

With the tight-binding band in the infinite- 
dimensional hypercubic lattice, divergence of the 
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Fig. 13. Temperature dependence of inverse staggered charge 
susceptibilities l/XQ Chg f° r conduction electrons exhibiting the 
CDW instability. The parameter J varies by 0.02 at every step 
to an adjacent curve. 



Fig. 15. Temperature dependence of inverse staggered charge sus- 
ceptibility l/Xq Chg around the quarter filling for J = 0.3. 
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Fig. 14. The CDW (q = Q) transition temperature Tcdw as a 
function of J for the quarter filling n c = 0.5. 



Fig. 16. The ground-state phase diagram of the infinite- 
dimensional Kondo lattice. The dashed lines indicate possible 
boundaries. 



Kondo temperature and the RKKY interaction. Hence, 
Doniach's picture is valid at half filling, although the 
estimation docs not take account of formation of the 
Kondo-insulating gap. We should note, however, that 
Doniach's picture fails in the regime far from half filling. 

The antifcrromagnetic state at half filling should be 
an insulator, since the Fermi surface with the nesting 
property disappears by the ordering. On the other hand, 
the paramagnetic regime at J > J c has been shown to 
be the Kondo insulator. Consequently, the ground state 
at half filling is an insulator regardless of the coupling 
constant J. The mechanism of the energy gap, however, 
differs between large- and small- J regimes. 



antifcrromagnetic susceptibility at Q = (tt, ■ ■ ■ , ir) has 
been demonstrated for a small-J case. With increasing 
coupling constant, the paramagnetic ground state takes 
over since antiferromagnctic fluctuations are suppressed 
by the Kondo effect. The critical coupling constant 
agrees with that estimated by comparison between the 



Ferromagnetism at low carrier densities 

We have confirmed divergence of the ferromagnetic 
susceptibility in the low-carrier-density regime. The ap- 
pearance of the ferromagnetism can be understood by the 
RKKY interaction. The variational method based on the 
mean-field theory has actually proposed its existence. 3 ' 4 
On the other hand, in the infinite-dimensional Anderson 
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lattice, Jarrell et al. did not find a divergent ferromag- 
netic susceptibility for any filling or set of parameters. 20 
Hence, the localized spins seem to be required for ferro- 
magnctism. 

The ferromagnetic phase extends beyond the critical 
coupling constant estimated from the comparison be- 
tween the Kondo temperature and the RKKY energy 
scale. This follows from the fact that all the local spins 
cannot be screened by a fewer number of conduction 
electrons. 8 The local susceptibility consequently follows 
the Curie- Weiss law even at the strong coupling. Our 
results imply that the effective interaction between the 
local spins is ferromagnetic. Note that the above argu- 
ment assumes the large number of interacting neigh- 
bors where the spin-pair correlation is less important. 
In low dimensions, on the other hand, intersite correla- 
tions among localized spins can make a collective singlet 
even without screening of conduction electrons as in the 
one-dimensional antiferromagnetic Heisenberg model. 

CDW state at quarter filling 

At quarter filling, we have found an instability of con- 
duction electrons against CDW formation at Q. The 
CDW ordering takes place in the middle range of the 
coupling constant. This work is the first confirmation of 
the CDW state in the Kondo lattice model to our knowl- 
edge. 

The origin of the CDW instability can be found in the 
strong-coupling effective model. In the strong-coupling 
limit, each conduction electron forms the Kondo singlet 
to screen the local spin. Since the spatial extent of the 
Kondo cloud is reduced to the lattice spacing at J = 
oo, the Kondo singlet may be regarded as a hole of the 
local spin. Accordingly, the strong-coupling limit of the 
Kondo lattice model can be mapped to the infinite-?/ 
Hubbard model, where the number of itinerant electrons 
is 1 — n c per site. In the effective Hamiltonian of order 
t 2 /J, we find a repulsion between the itinerant holes, 
or equivalently, between the Kondo singlets. 21,22 Hence, 
with decreasing J, the repulsion gets larger and may give 
rise to the CDW transition. This mechanism was first 
proposed by Hirsch for one dimension. 21 

Although the strong-coupling expansion predicts the 
CDW instability, whether the CDW ordering actually 
takes place or not depends on the range where the strong- 
coupling picture remains valid. Namely, if this picture 
breaks at J higher than the critical value for the CDW 
transition, the CDW state may not appear. Our results 
show that the strong-coupling picture holds away from 
J = oo in infinite dimensions. The appearance of the 
CDW state can be intuitively understood as the gain of 
the Kondo energy by the alternate distribution of con- 
duction electrons. 

Remaining issues 

In the present calculation, we have discussed instabili- 
ties only at the vector q = and Q. As has already been 
noted, there is a possibility of the ordering at the vec- 
tor other than q = and Q for an intermediate filling. 
In the mean-field theory at weak coupling, the ordering 
is expected at the peak of the polarization function H q , 



which can be incommensurate. It is an interesting issue 
whether the incommensurate ordering actually occurs or 
the commensurate vectors are chosen, when the local cor- 
relation is taken into account. 

The ground state in the CDW phase remains an open 
problem. Although occurrence of the CDW ordering has 
been demonstrated, the CDW state does not fix the spin 
configuration nor screen all the localized moments. In 
the strong-coupling picture, half the moments remain 
free at quarter filling. In order to address the mag- 
netic correlation in the ordered state, we need to work 
with a two-sublattice model. The problem is whether 
the ferromagnetic correlation observed outside the CDW 
phase remains below the ordering temperature. The two- 
sublattice calculations also provide informations on elec- 
tronic states in the antiferromagnetic phase as in the 
Hubbard model 23 and the Anderson lattice model. 24 Re- 
cently, the reduced antiferromagnetic moment in the 
Kondo lattice has been derived for the Bcthe lattice. 7 

Our formalism on the magnetic instabilities is based 
on the DMFT, so that the inter-site correlation is in- 
corporated in the mean-field level, which is justified in 
infinite spatial dimensions. In order to address two- or 
three-dimensional systems, spatial dependences of the 
self-energy as well as of the vertex part arc indispens- 
able. Extensions of the DMFT has been performed either 
by replacing the impurity with a cluster, 25 or by recon- 
structing the self-energy from the local vertex evaluated 
in the effective impurity. 26 ' 27 Investigation of the two- 
or three-dimensional Kondo lattice model is a most chal- 
lenging task, which should involve study of anisotropic 
superconductivity. 
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Appendix: The RKKY Interaction and Mean- 
Field Theory 

In this appendix, we derive the RKKY interaction, 
and estimate the magnetic transition temperature in the 
mean-field theory. Using the second-order perturbation 
theory for the Hamiltonian (1), we obtain the RKKY 
interaction of the form 

-Hrkky = - Jrkky(», j)Sj ■ Sj, (A-l) 

(ij) 

where the coupling constant Jrkky(«, j) is given by 

Jrkky(U) = 2J 2 n° j = 2 -fj2Ke iq<Ri - Rj) - (A-2) 

q 

The factor 2 comes from the spin degree of freedom. 

We analyze the RKKY model, eq. (A-l), in the mean- 
field theory. The spatial correlations of localized spins, 

(MF) . , 

Xq , are given by 

y ( mf ) = *L (A . 3 ) 

1 - JRKKYWXloc 
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where Jrkky(q') = 2J 2 n°. Without the renormalization 
effect by the conduction electrons, the susceptibility of 
the local spin follows the Curie law, X\ oc = 1/4T. A diver- 
gence of Xq MF ' ) yields the instability of the paramagnetic 
states. Consequently, the transition temperature Tq MF ' 
is determined by 



(A-4) 



We note that 11° in the right-hand side of this equation 
depends on the temperature. 
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